Observation of Topological Phase Transitions in Photonic Quasicrystals 
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Topological insulators and topological superconductors are distinguished by their bulk phase 
transitions and gapless states at a sharp boundary with the vacuum. Quasicrystals have recently 
been found to be topologically nontrivial. In quasicrystals, the bulk phase transitions occur in the 
same manner as standard topological materials, but their boundary phenomena are more subtle. In 
this work we directly observe bulk phase transitions, using photonic quasicrystals, by constructing 
a smooth boundary between topologically-distinct one-dimensional quasicrystals. Moreover, we 
use the same method to experimentally verify the topological equivalence between the Harper and 
Fibonacci quasicrystals. 



PACS numbers: 71.23.Ft, 73.43.-f, 73.43.Nq, 05.30.Rt 

The classification of gapped systems, such as band in- 
sulators and superconductors, by topological indices is a 
rapidly developing paradigm in condensed matter physics 
[TJ [2] . This novel approach provides insights to the char- 
acterization of states of matter, as well as predicts exotic 
phenomena. 

The topological classification of a system assigns an in- 
teger index to its energy gap. This index encodes prop- 
erties that are robust to distortions and deformation of 
the system. Hence, when a system with a given index is 
continuously deformed into a system whose index has a 
different value, the bulk energy gap must close, namely a 
quantum phase transition occurs. Accordingly, if a sys- 
tem can be continuously deformed into another system 
while keeping the bulk gap open, then their topologi- 
cal indices must be the same, defining them as topo- 
logically equivalent. Usually, such a phase transition 
manifests by the appearance of gap-traversing states at 
the boundary between topologically nontrivial materi- 
als and the topologically-trivial vacuum. Examples for 
such boundary states are the chiral modes of the inte- 
ger quantum Hall effect (IQHE), the Dirac cone of the 
three-dimensional topological insulator, and the Majo- 
rana fermions of one-dimensional (ID) topological super- 
conductors [TJ [2] . 

There is a variety of realizations to these exciting sys- 
tems in condensed matter, such as semiconducting het- 
erostructures |3j|4], Bismuth-based compounds [5HZ], and 
InSb nano- wires [8HTU]. Additional instances were ob- 
served in photonic crystals [HJH2], and suggested to ap- 
pear in cold atoms [15] . 

A new type of topological phenomena has been re- 
cently studied in quasiperiodic systems [HI [15]. Such 
systems, which are ordered but not periodic, were shown 
to be characterized by topological indices that are usu- 
ally attributed to systems of a dimension higher than 
their physical dimension. In particular, it was shown that 
the canonical ID quasiperiodic systems, i.e. the Harper 
(or Aubry- Andre) model [16] [17] , the diagonal Fibonacci 



model [18] [19], and their quasicrystalline off-diagonal 
variants [20, 21 , can be assigned Chern numbers. These 
number are topological indices that characterize generic 
two-dimensional (2D) systems [14] [15]. A continuous de- 
formation between two quasicrystals (QCs) with different 
Chern numbers will therefore result in a bulk phase tran- 
sition. Correspondingly, at a sharp boundary between 
such a QC and the vacuum, localized subgap boundary 
states may appear. These states were experimentally ob- 
served in photonic Harper QCs [14]. Further analyses 
of ID topological QCs with sharp boundaries have been 
conducted in cold atoms [22] [23] and superconducting 
wires [241J28]. 

The topological characterization of a QC is based on its 
long-range order, which is broken at the sharp boundary. 
Therefore, these boundary states do not always appear. 
However, if the boundary between the nontrivial QC and 
the trivial vacuum is formed by an adiabatic deformation, 
subgap states will always appear at the boundary, reveal- 
ing the bulk gap closure. Accordingly, such a deforma- 
tion can be used to prove equivalence between different 
quasiperiodic systems, if the energy gap remains open 
throughout the deformation. 

In this letter, we study the bulk gap closure that oc- 
curs when continuously deforming between topologically 
inequivalent quasiperiodic systems, and its absence when 
the systems are topologically equivalent. To this end, we 
create two inequivalent Harper QCs, and spatially in- 
terpolate between them. We then observe the closure of 
bulk energy gaps through the emergence of subgap states 
localized within the deformation region. In contrast, us- 
ing the same interpolation process between seemingly dif- 
ferent but topologically equivalent systems, the Harper 
and the Fibonacci QCs, no such phase transition occurs, 
thereby proving that these two models are indeed topo- 
logically equivalent. These phenomena are experimen- 
tally tested in quasiperiodic photonic lattices, where the 
phase transition, or its absence, are directly observed us- 
ing the propagation of light in waveguide arrays. 
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Photonic lattices are widely used for realizations of dif- 
ferent models originating from solid state physics, due 
to the high level of control over their parameters and 
behavior [29], [30]. Our photonic QCs are composed of 
an array of coupled single- mode waveguides, fabricated 
in bulk glass using femtosecond laser micro-fabrication 
technology [3TJ [32] . The overlap between the evanescent 
modes of the waveguides allows the propagating light to 
tunnel from each waveguide to its neighboring waveg- 
uides. Hence, the hopping amplitude between adjacent 
waveguides can be controlled by modulating the spacing 
between them. 

The dynamics of light propagating in these coupled 
waveguide arrays is described by the tight-binding model, 
with the propagation axis z taking over the role of time, 
id z i/j n = Hip n , where ip n is the wavefunction at waveg- 
uide number n. Taking the hopping amplitude to be real, 
we obtain the general Hamiltonian 



+ 1 > 



(1) 



where t n is the hopping amplitude from site n to site 
n — 1. 

Our intention is to study the transition that occurs 
when some system / is deformed into another system 
//, where each system has its own set of quasiperiodic 
hopping amplitudes t n and t n T , respectively. To this end, 
we fabricate a waveguide array with a deformed hopping 
profile t n = f n t n + (1 - f n ) tn, where 



1 

n—Li 
Ln 



1 < n < L/ 
Li < n < Li + L D , (2) 

Li + L D <n < Lj + L D + L n 



as depicted in Fig. [TJa). This procedure produces an 
array of length Li of system / on one side of the structure, 
an array of length Ln of system on the other side, 
and an L^-long deformation region, which continuously 
transforms between the two. This structure enables the 
study of the eigenstates of both systems as well as the 
transition between them, on a single waveguide array. 

The properties of the Hamiltonian fabricated within 
the photonic crystal are studied by injecting light into 
one of the waveguides in the array and measuring the 
outgoing intensity at the output facet using a CCD cam- 
era, as illustrated in Fig.Qb). The injected beam excites 
a wavepacket of all the modes that have a non-vanishing 
amplitude at the injection site, and the light propagates 
in the lattice according to this superposition of eigen- 
states. The width of the outgoing intensity distribution 
can therefore reveal the existence of localized eigenstates: 
If there is no localized state near the injection site, the 
light will spread freely throughout the array, propagating 
according to the bulk properties of the system. However, 
when light is injected in the vicinity of a localized state, 
its expansion is regulated by the width of the state. To 




(b) 



Coherent Light Source 



Lens 



m ► 



Lens 



D 



Waveguide 
Array 



(c) 



Ph 0.04 




FIG. 1. (Color online) Experimental methods, (a) Illustra- 
tion of a photonic waveguide array implementing a deforma- 
tion between two QCs [cf. Eq. Q]. (b) A schematics of the 
experimental setup. We focus a coherent light beam into a 
waveguide in the array, allow it to propagate along the struc- 
ture, and image the output intensity using a CCD camera, (c) 
Illustration of the relation between the localization parame- 
ter £ n and the measured intensity distribution \^ n \ 2 for two 
injection sites, n — 66 and n — 121. The shaded intervals de- 
note the A-distance neighborhood around the insertion point. 
When more light remains within this neighborhood, the value 
of £ n increases. 



quantify the localization of the outgoing light, we mea- 
sure the amount of light that remains within a small dis- 
tance A from the injection site n, by measuring the gen- 
eralized return probability [33] . 



in 



n+A Lj+Ln+Lj 

( E hM a )/( E 

m=n— A m—1 



l^nl 2 ) • (3) 



The relation between £ n and the intensity distribution 
\^ n \ 2 is illustrated in Fig. [ljc). Since £ n is meant to 
reveal the existence of localized states, we will choose A 
to be of the order of the width of a localized subgap state. 

Let us now introduce the specific quasiperiodic tight- 
binding models under study. In these models the hopping 
amplitude is modulated according to 



t n = t [1 + Xd n ] 



(4) 
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where to is the characteristic hopping amplitude of the 
system, A G [0, 1) is the modulation strength, and d n G 
[—1,1] is some quasiperiodic modulation function. 

Here we consider two such modulations: the Harper 
modulation 



dn = cos(27r6n + 
and the Fibonacci-like modulation 



(n + 2) 



(n + 1) 



1 = ±1. 



(5) 



(6) 



The long-range order of the Harper QC [20] is controlled 
by the modulation frequency b. Whenever b is irrational, 
the hopping modulation is incommensurate with the un- 
derlying lattice, resulting in a quasiperiodic pattern. The 
parameter <fi shifts the origin of the modulation. The en- 
ergy spectrum of the Harper model is composed of a frac- 
tal set of bands and gaps, in a way that depends on the 
modulation frequency b [20], [34]. Moreover, these gaps 
are associated with Chern numbers uniquely determined 
by b [15j [35]. Hence, two Harper QCs with bi ^ bjj 
are typically topologically inequivalent, since their gaps 
are associated with different Chern numbers. Compa- 
rably, the Fibonacci-like QC is constructed from a se- 
quence of two values that are ordered in a quasiperiodic 
manner. This sequence is obtained by applying the "cut- 
and-project" procedure on a square lattice onto the line 
y = x/r [36] . Whenever the slope of the line, r, is irra- 
tional, the sequence is quasiperiodic. For example, the 
case of r = (1 + Vb)/2 is the well-known Fibonacci QC 
[361137]. 

The properties of the Fibonacci-like QC differ in many 
ways from those of the Harper QC, e.g. the localization 
of the bulk wavefunctions [TBHUJ [38] [39]. Nevertheless, 
it was recently proven that they are topologically equiv- 
alent whenever the frequency of the Harper modulation 
satisfies b = (r + l)/r [15]. In such a case, the gaps 
of the Fibonacci-like QC are associated with the same 
Chern numbers as those of the Harper QC. Hence, for 
a given modulation frequency 6, the Harper QC can be 
continuously deformed into the Fibonacci-like QC with- 
out the appearance of a phase transition. 

Note that the deformation in Eq. ^ contains an ad- 
dition degree of freedom in the form of the parameter 
(j). This parameter has a crucial role in the observation 
of the topological boundary states of quasiperiodic sys- 
tems [13] , While the spectrum of our model is gapped 
in the bulk, localized boundary states appear, which tra- 
verse the energy gaps as a function of (j). Nevertheless, 
in this work we focus on bulk properties, which are </>- 
independent [40] . 

We now turn to our experimental results. We con- 
struct a deformation between topologically inequivalent 
Harper QCs with modulation frequencies bi ^ bn. Fig- 
ure |5|a) depicts the hopping amplitudes of a deformation 



between a Harper QC with bi = 2/(1 + \/5) and a Harper 
QC with b u = 2/(1 + \/CL5), where t 28/75 mm" 1 , 
A = 0.475 and <pi = <\>n = tt. The region lengths were 
taken to be Lj = Ljj = 84 and Ljj = 51. For this 
set of parameters, the bulk wavefunctions of the Harper 
QCs are extended [20 . We fabricated this structure in a 
75mm-long photonic waveguide array. This results in an 
effective propagation of 14 tunneling lengths, where the 
tunneling length is the characteristic length for hopping, 
namely 2 /to. With such a propagation length, light in- 
jected in the bulk of the structure will sufficiently expand 
in comparison to the width of localized subgap states. 

To experimentally observe the phase transition be- 
tween the two QCs, we injected light into each waveguide, 
and measured £ n as a function of the injection point n. 
The results are presented in Fig. |2|b). Two clear peaks 
in the deformation region can be seen over a relatively 
flat £ n outside the region. This is a clear indication of 
the existence of localized states within the deformation 
region. Note that measurements of £ n for n < 30 and 
n > 190 are omitted from this plot - for these injection 
sites the expanding light hits the edges of the structure, 
causing £ n to be skewed by boundary effects. 

To reveal the source of the peaks observed in £ n , we 
numerically obtain the local density of states (LDOS), 
which is presented in Fig. 2jc). The LDOS is defined 



by D n = (E - E m ) 



(m) 



where E m is the en- 

(m) 



ergy of the m th eigenstates, and ifn n) is its wavefunc- 
tion. D n (E) describes the spatial distribution of the 
eigenstates of the structure as a function of energy. For 
n < Li, we observe bands of extended states that cor- 
respond to the eigenstates of system /. Similarly, for 
n > Li + Ljj we recognize the band of extended states 
of system //. However, within the deformation region, 
there are few states that reside within the energy gaps 
and are spatially localized. These localized states tra- 
verse the energy gap along the deformation region, are 
an explicit signature of the phase transition between the 
inequivalent QCs [40], and are the origin of the measured 
peaks in £ n . We have therefore experimentally observed 
the bulk phase transition between two topologically in- 
equivalent Harper QCs. 

We now turn to study the transition between topo- 
logically equivalent QCs. We constructed a deformation 
between a Harper QC and Fibonacci QC, with a matched 
modulation frequency bi = (tjj + 1)/tjj = 2/(1 + y/S) 
[41]. The hopping amplitudes t n are depicted in Fig.|2|d), 
ior to = 28/75 mm" 1 , A = 0.225, (j) = tt(1 + 36), 
Li = Li i = 84 and Ljj — 51. For this set of parameters, 
the bulk wavefunctions of the Harper QC are extended, 
while those of the Fibonacci are critical [2T] , Neverthe- 
less, for the structure's effective propagation length, light 
injected into the bulk of both QCs will sufficiently expand 
in comparison to the width of potential localized states. 
The measured £ n of this system is depicted in Fig. [2|e), 
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FIG. 2. (Color online) Summary of results for a smooth deformation between (a-c) two topologically inequivalent Harper 
QCs, and (d-f) topologically equivalent Harper and Fibonacci QCs. In both experiments, Li — 84 (blue), Lu = 84 (red) 
and Ld = 51 (purple hues), (a) The hopping amplitude t n as a function of the lattice site n, for modulation frequencies 
bi — 2/(1 + y/b) and bu = 2/(1 + \/6.5) of the Harper QCs. (b) Experimentally- measured £ n for A = 7, as a function of the 
injection site n. The two peaks within the deformation region imply the existence of localized states, (c) Numerically obtained 
LDOS of the structure, D n (E). The energy bands are composed of extended states, while localized states (roughly 15 sites 
wide) traverse the gaps in the deformation region. These states manifest the transition between the inequivalent QCs. (d-f) 
Same as (a-c), but with a Harper QC deformed into a Fibonacci QC with bi = (tjj + 1)/tu = 2/(1 + Here, £ n shows no 
sign of localized states. Accordingly, though the distribution of the bands changes along the structure, no energy gap is closed. 
This proves the equivalence between the two QCs. 



showing no sign of localized states within the deforma- 
tion region. The numerically obtained LDOS is shown in 
Fig. [2^f ) . While the configuration of the bands changes 
considerably between the two QCs, no gap is closed along 
the deformation. Note, also, that two subgap states ap- 
pear at the sharp boundaries with the vacuum [40 . The 
open gaps and the corresponding absence of peaks in £ n 
serves as experimental evidence of the equivalence be- 
tween the Fibonacci and the Harper QCs. 

To conclude, in this paper we have presented a novel 
method to study topological phase transitions using a 
continuous deformation between two systems, which acts 
as a smooth boundary between them. When the bound- 
ary is sufficiently smooth, observations of subgap states 
localized within the deformation area serve as evidence 
of the phase transition. Such subgap states do not ap- 
pear when a phase transition does not take place, namely 
between topologically equivalent systems. Our method 
extends the prevailing approach which focuses on states 
that appear at sharp boundaries between topologically- 
nontrivial systems and the vacuum. It has proven use- 
ful in studying the topological characterization of ID 
quasiperiodic systems, since it circumvents the subtlety 
of their boundary phenomena at sharp boundaries. Fur- 



thermore, this technique may be useful to study other 
topological systems, such as (i) the weak and the crys- 
talline topological insulators, where the surface breaks 
the underlying symmetry [42H44] ; (ii) varying dopant 
concentration in 3D topological insulators [45 j; and (iii) 
nanowires that may host Major ana fermions at their 
boundaries [SHIP]. 
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SUPPLEMENTAL MATERIAL 



I. EFFECT OF FINITE-SIZED DEFORMATION 
ON THE PHASE TRANSITION 

In the main text, we study the phase transition that oc- 
curs when two topologically distinct quasicrystals (QCs) 
are deformed from one into the other, and its absence 
when the two QCs are topologically equivalent. This 
transition can be obtained by adiabatically deforming 
between the two QCs in real space, thereby creating 
a smooth boundary between them. As a result, when 
a topological phase transition occurs, the energy gap 
is continuously traversed by states localized within the 
deformation region. However, in physical systems, the 
deformation takes place over a finite length, making it 
not strictly adiabatic, and the gap closure is disrupted. 
Nevertheless, a finite number of subgap states remain, 
and their energies discretely traverse the gap. Here we 
show that these states serve as evidence of the topolog- 
ical phase transition, since their level spacing decreases 
with the deformation length. 

Figure [3] depicts the numerically obtained local density 
of states (LDOS) for structures deforming between two 
topologically distinct Harper QCs. For this plot, the pa- 
rameters are b T = 2/(1 + y/E), b H = 2/(1 + y/KE), Aj = 
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FIG. 3. Numerically obtained LDOS for deformation be- 
tween inequivalent Harper QCs, with bi = 2/(1 + y/E) and 
bn — 2/(1 + y/KE). (a) The deformation length is Ld — 60. 
The gap is traversed by 3 states in the deformation region, 
(b) For Ld — 300, there are 15 subgap states, and the level 
spacing between them decreases correspondingly. 



\jj 0.475, (\>i = (/>// = 0, and Lj + Ljj + L n 460. 
In Fig. |3ja), the deformation length is Ld = 60, and 3 
subgap states appear within each one of the large gaps. 
In Fig.|3jb), Ld = 300, and 15 subgap states appear. As 
Ld increases, the number of subgap states within a gap 
(^Vsub) increases, while the spacing between their energies 
(level spacing) decreases. Figure [4] plots the level spac- 
ing between two states residing in the middle of a large 
energy gap as a function of Ld- We can see that the 
level spacing scales as ~ 1/ Ld- Therefore, in the limit of 
Ld — >• 00, we approach the continuous phase transition. 




600 



FIG. 4. Level spacing between two states in the middle 
of a large gap as a function of the deformation length Ld- 
While the width of the states (l su b) is fixed, the number of 
states scales as Ld, and their level spacing scales as / su b /Ld • 
The dots are numerically obtained data, and the fitted curve 
matches 14.44/Ld- 

The level spacing can be approximated by A/7V su b, 
where A is the size of the gap. The scaling of 1/Ld 
implies that there is an additional length scale, namely 
the width of the subgap states, / su b. Given a fixed / su b, 
the number of subgap states within the deformation re- 
gion is approximately iV su b ~ Ld /hub, and their level 
spacing scales like Zsub/^D- The width / su b controls the 
smoothness of the boundary, and sets a lower bound on 
the deformation length for which subgap states always 
appear, i.e. Ld > hub- We can therefore conclude that 
structures with a finite deformation length larger than 
/ su b can indeed serve as a probe of the existence of a 
topological phase transition. 



II. BOUNDARY STATES AT A SHARP 
BOUNDARY AND THE ROLE OF 

We have seen that subgap states appear within a finite- 
length boundary between topologically inequivalent QCs. 
These states always appear for sufficiently long deforma- 
tion regions (Ld ^ ^sub), and their number increases 
with the deformation length Ld- 

In the opposite limit of Ld —> 0, the boundary be- 
comes sharp. At a sharp boundary, subgap states may 
or may not appear, depending on the exact pattern of the 
QCs [13] . Recall that the Harper modulation [cf. Eq. (5) 
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in the main text] includes the parameter 0. This param- 
eter shifts the origin of the quasiperiodic pattern, and is 
known as a phasonic degree of freedom. The appearance 
of subgap states at a sharp boundary with the vacuum 
is controlled by </>, since it determines the exact pattern 
of the QC at the termination point. As <j) varies from 
to 27r, boundary states continuously traverse the energy 
gaps. The number of gap-traversing states is given by 
the corresponding Chern number, which is determined 
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FIG. 5. Boundary states at a sharp boundary between two 
QCs. (a) Numerically obtained spectrum of a structure with 
the same parameters as in Fig. [3] and Ld = as a function 
of <j>i. The large gaps are traversed by two states when <j)i is 
scanned from to 2ty. (b) The corresponding LDOS for 0/ = 
0.17T. The gap-traversing state is localized at the boundary of 
QC I with the vacuum, (c) The corresponding LDOS for <j>i — 
IAtv. The subgap state is localized at the sharp boundary 
between QC I and QC II. (d) The spectrum of the structure 
as a function of Again, two states traverse the gaps, 

where one of them is localized at the boundary of QC II. The 
horizontal line in the middle of the gaps belongs the boundary 
state of QC /, which is unaffected by 



by the modulation frequency b. 

At a sharp boundary between two QCs, subgap states 
traverse the gap as either <pi or (\>n is scanned. This 
can be seen in Figs.|5|a)-(d), which show the numerically 
obtained spectra of a structure with the parameters from 
the previous section and Lb = 0, as a function of 
and (J)ji. As <pi is scanned, two states traverse each one 
of the large gaps, one localized at the sharp boundary 
of QC / with the vacuum, and the other at the sharp 
boundary between QC / and QC [as seen from the 
corresponding LDOS in Figs. [5jb) and|5|c)]. A similar 
behavior is observed when (j)u is scanned, only now, one 
of two the gap-traversing states is localized at the sharp 
boundary of QC // with the vacuum. 

In contrast to a sharp boundary, for infinitely long Ld , 
the gap traversing states are independent of <\>. For fi- 
nite Ld, where there is a finite number of subgap states, 
the gap completely closes when <fi is scanned. This can 
be seen in Fig. [6j which depicts a numerically obtained 
spectrum with Ljj = 80 as a function of </>/. As <\>i is 
scanned from to 27r, the energy of each subgap state 
in the deformation region continuously connects to the 
energy of the next level. We can also recognize the state 
at the boundary of QC / with the vacuum [cf. Fig. [5ja)]. 
Notably, the motion of the eigenstates with <\>i represents 
pumping analogues to Laughlin's pumping in the integer 
quantum Hall effect 
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FIG. 6. Spectrum of a structure with the same parameters 
as in Fig. [3] and Ld — 80 as a function of <fii. The energy of 
each subgap state in the deformation region approaches the 
energy of the next one when <pi is scanned. The state at the 
boundary of QC / with the vacuum behaves like in Fig.[5ja). 



